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(Dated: February 7, 2008) 

The first observation of an azimuthal cross-section asymmetry with respect to the charge of 
the incoming lepton beam is reported from a study of hard exclusive electroproduction of real 
photons. The data have been accumulated by the HERMES experiment at DESY, in which the 
HERA 27.6 GeV electron or positron beam scattered off an unpolarized hydrogen gas target. The 
observed asymmetry is attributed to the interference between the Bethe-Heitler process and the 
Deeply Virtual Compton Scattering (DVCS) process. The interference term is sensitive to DVCS 
amplitudes, which provide the most direct access to Generalized Parton Distributions. 

PACS numbers: 13.60.-r, 24.85. +p, 13.60.Fz, 14.20.Dh 



The partonic structure of the nucleon has been tradi- 
tionally described in terms of Parton Distribution Func- 
tions (PDFs), which appear in the interpretation of, 
e.g., inclusive Deeply Inelastic Scattering (DIS). More 
recently, PDFs have been subsumed within Generalized 
Parton Distributions (GPDs) 0, 0, 0, which relate di- 
rectly to hard exclusive processes that involve at least 
one additional hard vertex, yet leave the target nucleon 
intact. The ordinary PDFs and nucleon elastic form fac- 
tors appear as kinematic limits and moments of GPDs, 
respectively. Strong interest in the formalism of GPDs 
has emerged after GPDs were found to offer the first pos- 
sibility to reveal the total angular momentum carried by 
the quarks in the nucleon [2j. More recent discussions 
focus on the potential of GPDs as a three-dimensional 
representation of hadrons at the parton level 0, 0, @, 0, S] - 

Among all practical probes, the DVCS process, i.e., 
the hard exclusive leptoproduction of a real photon 
(ep — > epj), appears to provide the theoretically clean- 
est access to GPDs. Direct access to the DVCS ampli- 
tudes is provided by the interference between the DVCS 
and Bethe-Heitler (BH) processes, in which the photon 
is radiated from a parton and from the lepton, respec- 
tively. Since these processes have an identical final state, 
the squared photon production amplitude is given by 

M = |T Bff | + |t dvcs | + T DVCS T* H + Tp VCS T BHl (1) 



where I denotes the interference term. It introduces a 
dependence on the beam charge, which is a rare phe- 
nomenon normally confined to processes involving the 
weak interaction. The BH amplitude (t bh ) is precisely 
calculable from measured elastic form factors. The cross 
section depends on the Bjorken scaling variable xb, the 
squared virtual-photon four-momentum —Q 2 , and the 
squared four-momentum transfer t to the target. This 
interpretation of the virtual-photon kinematics does not 
apply to the BH process. In addition, the cross section 
depends on the azimuthal angle <fi S [— 7T, 7t], denned as 
the angle between the plane containing the incoming and 



outgoing lepton trajectories and the plane correspond- 
ingly defined by the virtual and real photon 

For an unpolarized proton target, and —t « Q 2 , the 
interference term is given by [lfj| 

I (x -C[a cos^ReTW 1 ' 1 -bP t sin^ImM 1 ' 1 ], (2) 

where the lepton beam has longitudinal polarization P; 
and charge C = ±1, and a and b are functions of the ra- 
tio of longitudinal to transverse virtual-photon flux. A 
polarization-independent constant term, higher harmon- 
ics (cos 2(f>, cos3</>, sin 20), as well as a cos</> dependence 
in the prefactor, have been neglected since they are sup- 
pressed by at least 0(l/Q) or 0(a s ). Here a s is the 
strong coupling constant. The squared BH and DVCS 
amplitudes have their own <j) dependences, but do not 
depend on the sign of the charge. Hence the measure- 
ment of a cross section asymmetry with respect to the 
beam charge (da + — da~) / '(da + + da~) is a way to sin- 
gle out the interference term [llj in the numerator, while 
the denominator is dominated by a ^-independent BH 
contribution. See Ref. [l2| for details and full equations. 
The photon-helicity-conserving amplitude 
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is given by a linear combination of the Compton Form 
Factors (CFFs) H, H and £, to geth er with the Dirac 
and Pauli form factors Fi and F^ 12[ . Here M p denotes 
the proton mass. The CFFs are convolutions of the cor- 
responding twist-2 GPDs H, H and E with the hard 
scattering amplitude. 

The sin <p modulation accessing the imaginary part of 

M. 1 ' 1 has already been observed in cross-section asym- 
metries with respect to the beam helicity ljl, 14 j. This 
Letter reports the first measurement of an asymmetry 
with_respect to the beam charge, accessing the real part 
of M. 1,1 via a cos</> modulation. 

Data with an unpolarized hydrogen target were accu- 
mulated using the HERMES spectrometer and the 
longitudinally polarized 27.6 GeV electron and positron 
beams of the HERA accelerator at DESY. Events were 
selected if they contained exactly one photon and one 
charged track identified as the scattered lepton. The 
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FIG. 1: Distributions in missing-mass squared from data (sta- 
tistical error bars) and from Monte Carlo simulations (line). 
The latter include elastic BH (filled area) and associated BH 
(hatched area) processes as well as semi-inclusive background. 
The vertical lines enclose the selected exclusive region. 



FIG. 2: Beam-charge asymmetry Ac for the hard electropro- 
duction of photons off protons as a function of the azimuthal 
angle \<j>\, for the exclusive sample before background cor- 
rection. Statistical uncertainties are shown. The solid curve 
represents the four-parameter fit: (-0.011 ±0.019) + (0.060± 
0.027) cos <f> + (0.016 ± 0.026) cos 20 + (0.034 ± 0.027) cos 30. 
The dashed line shows the pure cos dependence. 



hadron contamination in this lepton sample was kept be- 
low 1% by combining the information from a transition- 
radiation detector, a preshower scintillator detector, and 
an electromagnetic calorimeter. The kinematic require- 



ments imposed were 1 GeV < 



< 



10 GeV 2 , 0.03 



<x B < 0.35, W > 3 GeV, and v < 22 GeV, where W de- 
notes the initial photon-nucleon invariant mass and v is 
the virtual-photon energy in the target rest frame. The 
real photon was identified by detecting an energy depo- 
sition above 5 GeV in the calorimeter in addition to a 
signal in the preshower detector, without an associated 
charged track in the back region of the spectrometer. Un- 
like in a preliminary stage of this analysis 16], the useful 
range of the polar angle (9 7 » 7 between the virtual and real 
photons is not limited by the photon position resolution, 
due to an improved reconstruction algorithm. Now the 
restriction on # 7 * 7 is relaxed to # 7 * 7 > 5 mrad, limited 
mainly by the electron momentum resolution. In addi- 
tion, a stricter upper limit of # 7 * 7 < 45 mrad is imposed 
in order to improve the signal-to-background ratio 17 1 . 

The recoiling proton was not detected. Hence exclusive 
events are selected by requiring the missing mass Mx of 
the reaction ep — > cyX to be close to the proton mass. 
Fig.[T]shows the distribution in M x = {q+p—q') 2 , with q, 
p, and q' being the four-momenta of the virtual photon, 
the target nucleon in the initial state, and the real pho- 
ton, respectively. Mainly due to the resolution in photon 
energy, the exclusive peak extends to negative values, in 
which case Mx is defined as — y/— M x . The exclusive 
region is defined as (-1.5 GeV) 2 < M x < (1.7 GeV) 2 , 
based on the result of a Monte Carlo simulation (MC), 
shown in the same figure. The Mo and Tsai formalism 
2l| is used to simulate the elastic BH process leaving 



the target nucleon intact, and the associated BH pro- 
cess, where the nucleon is excited to a resonant state. 
For the latter a cross-section parametrization for the res- 



onance region is used [22j. Not included in the simula- 
tion is the DVCS contribution, which in this kinematic 
regime is ex pec ted to be much smaller than that of the 
BH process [23j. The simulation also takes into account 
the semi-inclusive production of neutral mesons (mostly 
7T°), where all but one of the decay photons escape de- 
tection. For this, the MC generator LEPTO in con- 
junction with a special JETSET [19] fragmentation tune 
is used, the latter being optimized for energies relevant to 
HERMES [13]. Not shown in Fig. Q] is the contribution of 
exclusive tt° production, which contributes less than 2.5% 
to the exclusive region based on the model in Ref. [2||. 
As shown in Fig. [TJ data and MC are in good agree- 
ment taking into account that they are both absolutely 
normalized, and that the MC does not include radiative 
corrections to the BH cross-section. The MC yield ex- 
ceeds the data by about 20% in the exclusive region, as 
may be expected (25| . A simulation including second or- 
der radiative processes should give an improved MC-data 
comparison in the full missing-mass region, since a part 
of the exclusive events experiencing second order radia- 
tive corrections will not be lost but reconstructed in the 
non-exclusive region. 

The beam-charge asymmetry is evaluated as 



A c {(j>) = 



JV+(0)-JV-(0) 
N+(<f>)+N-{<f>y 



(4) 



where N + (cf)) and N~(cf>) represent the single-photon 
yields per bin, normalized to the number of detected in- 
clusive DIS events using the positron and electron beam, 
respectively. Since these beams were polarized, sinu- 
soidal contributions appear in numerator and denomina- 
tor of Eq.[4j with the last term in Eq.[2]giving the biggest 
contribution. In order to cancel these contributions, the 
'symmetrized' beam-charge asymmetry is calculated by 
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FIG. 3: The cos (j> amplitude of the beam-charge asymmetry 
as a function of the missing mass, before background correc- 
tion. Statistical uncertainties are shown. 

replacing (f> by in Eq. [U The result for the exclusive 
region is displayed in Fig. [5] The shown four-parameter 
fit yields a non-zero cos (j> amplitude of 0.060±0.027. The 
result after the background correction described below is 
given in Table HI last row. The constant term as well 
as the cos 20 and cos 30 terms are compatible with zero. 
Fig. [3] shows the cos</f> amplitude in several Mx bins. At 
higher Mx the result is compatible with zero, confirming 
the absence of charge-dependent instrumental effects. 

As the recoiling proton remains undetected, t is in- 
ferred from the measurement of the other final-state par- 
ticles. For elastic events, kinematics relate the energy 
with the direction of the real photon, opening the possi- 
bility to omit the real-photon energy, which is the quan- 
tity subject to larger uncertainty. Thus the value of t in 
the exclusive region is calculated as 



2 v [v — \Jv 2 + Q 2 cos 9 1 * 1 ) 



1 + ir ( v ~ V" 2 + Q 2 cosG 



(5) 



The error caused by applying this expression to inelastic 
events (~17% in the exclusive region) is accounted for in 
the MC simulation that is used to calculate the fractional 
contribution of background processes per kinematic bin 
in — t (see Ref. 17] for details). 

Figure [5] shows the cos <j) amplitude derived from the 
four-parameter fit in each of four bins in — t. In each bin, 
this result is corrected for the semi-inclusive background, 
which is treated as a dilution since the background asym- 
metry can only be non-zero at next-to-leading order in 
QED. The total background contribution is about 6% as 
derived from the MC simulation, wherein the elastic and 
associated BH contributions are scaled down by the 20% 
described above. The resulting cos0 amplitudes are ex- 
pected to originate from only elastic and associated pro- 
duction. The associated BH processes contribute about 
5, 11, 18, and 29% to the yields in the four —t bins, or 
11% in the full t-range, with an estimated fractional un- 
certainty of 10%. The dominant contribution to the to- 
tal systematic uncertainty of the cos (f> amplitudes stems 
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FIG. 4: The cos 4> amplitude of the beam-charge asymmetry 
as a function of —t for the exclusive region (-1.5 GeV < Mx < 
1.7 GeV), after background correction. The error bars (band) 
represent(s) the statistical (systematic) uncertainties. The 
calculations based on GPD models 26, 27] use either a factor- 
ized i-dependence with (dashed-dotted) or without (dotted) 
the D-term contribution, or a Regge-inspired t-dependence 
with (dashed) or without (solid) the D-term contribution. 



from effects due to possible deviations of the detector 
and/or the beam from their nominal positions. These 
effects can be as large as 0.02 per bin. Based on the 
models in Ref. 23[ , acceptance and smearing effects can 
contribute up to 20% of the cos0 amplitude, and thus 
dominate the systematic uncertainty in the last —t bin. 
The other sources of systematic uncertainties are due to 
a possible difference in the calorimeter calibration be- 
tween the two data sets, the uncertainties from the semi- 
inclusive background correction described above, and the 
dilution of the asymmetry due to exclusively produced 
7T° mesons misidentified as photons. These contributions 
are combined quadratically in the total systematic uncer- 
tainty per bin in —t, given in Table [If. Not included is 
any contribution due to additional QED vertices, as the 
most significant of these has been estimated to be negligi- 
ble, at least in the case of polarization asymmetries 24 1. 

The theoretical calculations for the ep — > epj pro- 
cess shown in Fig. U employ GPD models developed in 
Refs. 26[ 27 1, which are based on the widely used frame- 
work of double distributions 28] . The model parameters 
of interest are those that change the GPD H since the 
impact of the GPDs H and E is suppressed at small val- 
ues of xb and —t, respectively (cf. Eq. [3]). The code of 
Ref. [2911 was used to calculate the values for the cos0 



1 Note that a preliminary result of this analysis [l6| with a t— 
averaged value of 0.11 ± 0.04 (stat.) ±0.03 (sys.) was derived 
at a much larger mean —t value ({— t) = 0.27 GeV 2 ) due to 
different requirements on 7 » 7 , as described above. It thus can- 
not be compared to the t-averaged result given in Table [I] but 
approximately to the result in the third — t bin. 
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—t bin 


i-t) 


(x B ) 


(Q 2 ) 


A™ s0 ± stat. ± sys. 


(GeV 2 ) 


(GeV 2 ) 




(GeV 2 ) 




< 0.06 


0.03 


0.08 


2.0 


0.024 ± 0.043 ± 0.022 


0.06 - 0.14 


0.09 


0.10 


2.6 


0.020 ± 0.054 ± 0.022 


0.14 - 0.30 


0.20 


0.11 


3.0 


0.071 ± 0.066 ± 0.028 


0.30 - 0.70 


0.42 


0.12 


3.7 


0.377 ± 0.110 ± 0.081 


< 0.70 


0.12 


0.10 


2.5 


0.063 ± 0.029 ± 0.028 



TABLE I: The coscj) amplitude of the beam-charge asymme- 
try per kinematic bin in —t after background correction and 
the respective average kinematic values. 



the first time. A cos <j> dependence has been observed in 
a kinematic region where the target proton is predom- 
inantly left intact. The cos0 dependence is predicted 
to arise from the interference between the deeply virtual 
Compton scattering and Bethe-Heitler processes. These 
data can already be used to distinguish among theoretical 
models for generalized parton distributions. 

We gratefully acknowledge the DESY management for 
its support, the staff at DESY and the collaborating insti- 
tutions for their significant effort, and our national fund- 
ing agencies for financial support. 



amplitude of the beam-charge asymmetry at the average 
kinematics (see Table H]) of every —t bin and not at the 
kinematics of every event since it is too computationally 
intensive. The difference between these two approaches 
is strongly model dependent: Tests [30j show differences 
of up to 20% using the models in Ref. [23j, which are 
equivalent to the factorized models in Fig. 2J Four dif- 
ferent parameter sets are selected by choosing either a 
factorized or a Regge-inspired i-dependence, each with 
or without the contribution of a negative value of the so- 
called D-term [31], which is related to the spontaneous 
breaking of chiral symmetry in QCD 27]. It contributes 
to the real part of the DVCS amplitude only and there- 
fore can be investigated with the beam-charge asymme- 
try for the first time. A large positive D-term is ruled 
out by these data [16j, since it leads to a negative value 
for the beam-charge asymmetry. The parameters b va i 
and b sea in the profile function are fixed at unity, 
since the beam-charge asymmetry has been shown to be 
largely insensitive to them [l7|. In comparing the pre- 
dictions to the data at large —t, it should be borne in 
mind that the model calculations do not include associ- 
ated production, which increases with —t as mentioned 
above. The three data points at small —t exclude the 
model based on the Regge-inspired i-dependence with 
the D-term contribution. 

These data are in agreement with a very recent cal- 
culation based on a dual parametrization of GPDs and 
either of two models for the i-dependence [33[. Like the 
calculations shown in Fig. [4] these calculations were per- 
formed at the average kinematics (see TableHJ of every —t 
bin. Earlier model calculations were mostly done not as a 
function of t but at the average kinematics of the prelimi- 
nary result [l6| , and thus cannot rigorously be compared 
to this measurement, see, e.g., Refs. [7J, |L2|, |3J, |35| . How- 
ever, since they span a wider or different range for the 
magnitude of the beam-charge asymmetry when com- 
pared to the models shown here, it is apparent that this 
measurement can constrain the GPD H . 

In conclusion, a beam-charge azimuthal asymmetry in 
electroproduction of real photons has been measured for 
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